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INTRODUCTION 

The  American  Association  of  State  Highv/ay  Officials' 
PolJ-C'/  on  Geometric  Design  states  that  "the  appearance  aspect 
of  superelevation  runoff  largely  governs  its  length"  (1) 
when  referring  to  spirals  (transition  curves)  and  their  use. 
Likev/ise,  Pushkarev  (2)  proposes  "that  a  frankly  esthetic 
approach  to  transition  curves  is  justified,  with  their  length 
determined  not  by  psuedo-utilitarian  minima,  but  rather  by 
what  is  visually  necessary  to  achieve  a  generous,  free-flowing 
continuity  of  alignment." 

Cron  (3)  adds  to  this  by  an  analogy  to  the  railroad 
designer. 

The  railroad  locators  liked  to  alternate  right  and 
left  curves,  but  were  also  careful  to  keep  a  good 
length  of  tangent  between  them.   This  left  ample 
space  for  the  later  insertion  of  spiral  transitions 
into  the  alignment.   Highway  designers,  likewise, 
should  provide  ample  tangents  betv;een  curves  in 
opposite  directions,  and  should  also  provide  for 
spirals  at  the  beginning  and  ending  of  all  but  the 
very  flattest  horizontal  curves.   The  spiral  helps 
the  driver  to  stay  in  his  lane  when  entering  or 
leaving  a  curve;  it  provides  a  convenient  and  mathe- 
matically correct  way  to  superelevate  or  "bank"  the 
curve;  and  it  also  greatly  improves  the  appearance 
of  the  highway,  particularly  where  the  edges  are 
sharply  defined,  as  in  concrete  pavements.   Since 
the  application  of  spirals  costs  nothing  except  a 
little  more  figuring  during  the  location  survey, 
their  use  can  be  fully  justified  on  aesthetic 
grounds  alone. 

Indeed,  the  addition  of  spirals  to  the  horizontal  align- 
ment of  a  roadway  provides  a  measurable  visual  benefit  at  a 
very  small  increase  in  cost.   The  main  endeavor  of  this 
research  was  to  develop  criteria  for  the  selection  of  spiral 


lengths  based  on  the  visual  appearance  of  the  curve  so  this 
benefit  can  be  realized. 


LITERATURE  SEARCH 

In  addition  to  the  previously  mentioned  references,  the 
problem  of  highway  alignment  coordination  has  been  recognized 
by  numerous  other  authors  (4,  5,  6).   Most  of  the  discussion 
of  the  problem,  however,  was  limited  to  general  statements 
concerning  the  visual  appearance  of  the  roadway.   Although 
the  comments  were  general  in  nature,  they  did  point  out  areas 
of  possible  problems  in  highway  alignment. 

Smith  and  Yotter  (7)  utilized  perspective  drawings  to 
determine  the  minimum  acceptable  length  of  sag  vertical  curves 
that  would  provide  an  aesthetically  pleasing  roadway.   They 
also  investigated  the  problem  of  a  small  change  of  direction 
in  the  horizontal  alignment.   Their  research  showed  that  the 
distance  from  which  a  curve  is  viewed  and  the  angle  from  which 
it  is  viewed  affects  the  appearance  of  the  curve.   Because  the 
main  concern  of  these  researchers  was  the  appearance  of  a 
sag  vertical  curve,  the  results  of  their  research  were  of 
limited  value  in  the  detailed  search  for  selection  criteria 
of  spiral  curves. 

T.  Ten  Brumjnelaar  (8)  advocates  the  use  of  deflection  and 
curvature  diagrams  to  determine  locations  which  can  cause 
discomfort  or  hazard  in  roadway  design.   The  technique  involves 
constructing  the  equivalent  to  shear  and  moment  diagrams  for 
the  deflection  and  curvature  of  tlie  roadway  alignment.   The 
shape  of  these  diagrams  will  show  areas  which  are  not  geometrically 
consistent  with  the  adjoining  sections  of  roadway.   The  use  of 
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these   diagrams    seems   to  be    limited  to  the   role   of   comparison 
of   sections   of   one  .highv/ay   as   no   indication   of   the  maximum 
allov;able  param.eters  were    given, 

Godin  et   al    (5)    stated  that   the    length   of  spiral   cuirve 
required  for  visual   appearance   exceeded  both   the   length  needed 
for  superelevation   runoff   and  the    length  needed  to   allow   a 
driver  to  steer  a   smooth   transition   from  tangent  to  circular 
curve , 

Only  two   references    contained  quantitative   inform.ation 
concerning  desirable    lengths   of   spiral   curves  based  on   the 
visual  aspect   of  the    roadv;ay.      Pushkarev    (2)    recommends   that 
the   spiral   length   should  have   a   ratio  to  the   circular  curve 
of   1:2:1,      He   further  states   that  spirals  whose    length  becomes 
too  great  in   relation    to  the    circular  curve   length   cause    the 
total   curve   to  appear   to  have   a   sharp  bend  in   the  middle. 
Godin    (9)    recommended  that  the  minimum  length   of   spiral   curve 
should  be   R/9 ,  v/here   R  is    the   limiting  radius   of   the  spiral. 

In   a  search   for   additional   specific,   quantitative   criteria 
concerning  v/hat   constitutes    a  visually-pleasing  highway   align- 
ment,  psychology  writings  were   investigated  for  a  clue   as   to 
what   the  human   eye    "sees"    as   pleasing  and   for  what  reason, 
VJritings    on   Gestalt  Psychology,    in  particular,   v/ere   examined 
(10,11).      These   investigations   failed  to  uncover   any  useful 
quantitative   results. 


PURPOSE 

The  purpose  of'  the  research  described  in  this  paper  was 
to  determine  criteria  for  the  selection  of  a  spiral  length 
for  horizontal  curves  to  aid  the  highway  designer  in  the 
creation  of  a  visually  pleasing  roadway. 


SCOPE 


The   research  was    limited  to: 

1,  The  horizontal    curve  v/as   studied,   with   and  v/ithout  a 
spiral   curve,    to  determine   the    relationship  between   the   sighting 
distance,    the   geometry   of    the   curve   and  the  angle   from  which 
the   curve   is   viewed. 

2,  The   length   of   spiral   curve  necessary  to  provide   a 
pleasing   appearance  was    investigated. 


METHOD  OF  SOLUTIQN 

From  personal  observations   and   the   information   gained   from 
the   literature   search,    it  was  hypothesized  that  the    factors 
affecting  the   visual  appearance   of   a  horizontal   curve  were 
sight  distance    (SD)    display   angle    (DA)    and  the   geometry   of   the 
curve.      Sight   distance   is    defined  as    the   distance    from  ttie 
observer  to  the  beginning  of   the   curve,   PC  or  TS  v/hichever  is 
appropriate;    display   angle   is   defined  as    the   angle  between  the 
observer's   line   of   sight  and  the  plane   containing  the  PC   or  TS 
of   the   curve;    and  the   geometry   of   the   curve   is   defined  as 
the   spiral   length,    if   any,    and  the    limiting  radius   of   the   curve. 
Figure    1  gives    a  pictorial  definition   of   sight  distance   and 
display   angle. 

Due   to  the   almost  total    lack   of  existing  highways   contain- 
ing spiral   curves   as    a  design  element,    it  was    felt   that  various 
geometric   conditions  would  have  to  be   simulated  to  provide 
sufficient  data   for  analysis   of   the   problem  of   selecting  spiral 
curve   lengths    for   a  pleasing  visual  appearance. 

A  computer  algorithm    (7)    was   available   that  v;ould  convert 
three-dimensional   coordinates    into   two-dimensional   coordinates 
which,   when   dravm,    gave    a  perspective  view.      Figure   2    gives   a 
graphic  illustration   of   this    algorithm.      The  three-dimensional 
coordinates    of   the   roadway  were   used  as    inputs.      Then,    after 
defining  the   center  of   interest   coordinates   and   the   observer 
position   coordinates,   both   in  three   dim.ensions,   the   tv;o- 
dimensional  perspective   coordinates  were   determined  by   tlie 
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projection   of  the   roadway   coordinates   on  a  plane  perpendicular 
to  the    line   of   sight   from  the   observer  to   the   center  of   interest, 
The  projections   on  this  plane,    called  the  picture  plane,    form  a 
perspective   view  of   the   three-dimensional   roadway   coordinates. 
The  equations    for   calculating  the   perspective  plane   coordinates 
are   included   in   Figure   2.      The   ready  access   to  this   algorithm 
was   instrumental   in   the   selection   of   this   approach   to  the 
problem  of   selection   of  spiral  curve   lengths   for  visual  appear- 
ance. 

The   only   remaining  obstacle   to  the   utilization   of  this 
approach  was   the   generation  of  the   three-dimensional   roadway 
coordinates.      This  was   more   of  a  problem  than  was   originally 
anticipated.      It  was  hoped   that  the   coordinate   GeOm.etry    (COGO) 
portion   of  the   Integrated  Civil  Engineering  System    (ICES) 
developed  by  the   Civil   Engineering  Systems   Laboratory   at  the 
Massachusetts   Institute   of   Technology    (MIT)    could  be   used  to 
generate   the   three-dimensional   coordinates.      The   ICES   system 
(12)    available    for  the   IBM   360    computer  included  provision 
for  the   calculation   of  spiral   curves.      However,   after  discussion 
with   computing  personnel   at  Kansas   State  University   and   the 
persons   in   charge   of   ICES    at  MIT,    it  was   found  to  be  extremely 
difficult   to  secure   any   form  of   output  other  than   the   standard 
printed  output.      This  v/ould  have  necessitated  punching   all  of 
the   three-dimensional   coordinates   for  input   into  the   coordinate 
transformation   program. 

Therefore,    it  was   deemed  necessary'  to  develop   a   computer 
program  to  generate   the    roadway   coordinates.      The  equations 
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for  calculating  the  centerline  roadway  coordinates  are  froni 
Hickerson's  book  (13)  Route  Surveys  and  Design,   The  spiral 
curve  used  was  approximately  a  cubic  parabola.   The  general 
equation  for  the  spiral  is  R  =  K/1 ;  where,  R  is  the  radius  of 
the  spiral  at  any  point  on  the  spiral,  1  is  the  distance  along 
the  spiral  to  the  point  of  radius  R  and  K  is  a  constant.   From 
this  equation  it  can  be  sho\m  that  the  radius  of  the  spiral  is 
infinite  v/hen  1  =  0,  at  the  TS  or  beginning  of  the  spiral  curve, 
and  decreases  as  1  increases. 

The  equations  used  for  calculation  of  the  spiral  curve  are 
shown  in  Figure  3,   The  equations  for  x  and  y   are  very  close 
approximations  of  the  true  equations,  and  the  loss  of  accuracy 
due  to  these  approximations  was  not  felt  to  be  sufficient  to 
alter  the  perspective  picture. 

The  edge  of  roadway  and  edge  of  shoulder  coordinates  were 
calculated  by  determining  the  centerline  direction  and  using 
sine  and  cosine  functions  to  locate  them  at  a  given  offset 
distance.   Again,  this  approximation  to  the  true  coordinates 
of  these  points  was  not  felt  to  be  great  enough  to  alter  the 
realism  of  the  perspectives. 

The  resulting  program  allows  any  roadv;ay  geometries  to  be 
sim.ulated  with  a  minimum  of  input  data.   It  can  be  used  to 
simulate  a  hypothetical  situation  or  give  the  coordinates  of 
an  actual  location.   The  coordinate  generation  program  was 
then  made  a  subroutine  of  the  coordinate  transformation  and 
plotting  program.   The  program  output  was  a  printout  of  the 
centerline  geometry  and  data  stored  on  a  7-track  comouter 
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FIGUIIE  3.    Spiral  Curve  Calculation  Equat5.ons 
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tape  which  could  be  mounted  on  a  Calcorap  Incremental  Plotter 
for  drawing  the  perspective  pictures.   A  general  flow  diagrara 
of  the  program  is  presented  in  Figure  4.   A  printout  of  the 
complete  program  is  included  in  the  appendix. 

Simulation  of  roadway  geometry  would  be  valid  only  if 
the  perspective  picture  looks  like  the  actual  location.   The 
realism  of  the  perspective  pictures  can  be  judged  by  comparing 
the  perspective  drawing  and  photograph  of  the  same  location  in 
Figure  5.   The  observer  position  for  the  perspective  picture 
was  approximately  the  same  location  from  v;hich  the  photograph 
was  taken.   The  photograph  was  taken  with  a  35mm,  NIKON  F, 
automatic  single  lens  reflex  camera  with  a  NIKKOR,  f3.5,  43mm 
to  86mra  Zoom  lens.   The  lens  was  set  at  86mm,  the  focal  length 
judged  to  provide  the  most  "natural"  perspective.   It  should 
be  noted  that  only  the  right  hand  lanes  of  the  divided  highway 
are  represented  in  the  perspective  drawing,  since  there  v/as  no 
provision  for  drawing  a  divided  highway  in  the  original  program. 

Approximately  five-hundred  perspective  pictures  were 
plotted  with  varying  sight  distance,  limiting  degree  of  curve, 
spiral  length  and  display  angle.   For  any  given  set  of  geometric 
conditions,  the  sight  distance  and  display  angle  v;ere  varied 
by  moving  the  observer  position.   Moving  the  observer  up  into 
space  above  the  roadway  created  artificial  display  angles  and 
simulated  the  case  v;here  a  sag  vertical  curve  is  located 
between  the  observer  and  the  horizontal  curve  without  the  dis- 
ruptive element  caused  by  the  vertical  curve.   It  was  felt  that 
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FIGURE  5.   Photograph  and  Perspective  of  Selected 
Location 
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if  a  vertical  curve  had  been  used  to  create  a  display  angle, 
the  appearance  of  the  horizontal  curve  might  be  affected  by 
the  vertical  curve,  thus  creating  a  problem  of  alignment 
coordination. 

A  wide  range  of  geometric  conditions  was  investigated. 
The  degree  of  curve  was  varied  from  30  minutes  to  5  degrees. 
Some  of  the  curves  were  constructed  without  a  spiral  transi- 
tion so  that  the  visual  effect  of  adding  a  spiral  curve 
could  be  tested.   The  spiral  lengths  varied  from  100  to  2000 
feet.   An  attempt  was  made  to  select  lengths  which  would  test 
the  theories  proposed  by  Godin  and  Pushkarev.   The  lengths 
of  spirals  examined  at  each  limiting  degree  of  curve  are 
shown  in  Table  1, 

All  curves  were  rated  according  to  their  smoothness  of 
appearance  and  the  rate  at  which  they  seemed  to  diverge  from 
the  tangent.   The  rating  scale  was:   acceptable,  questionable 
or  unacceptable.   Although  this  may  seem  a  crude  rating  scale, 
it  was  felt  that  any  refinement  beyond  this  level  was  not 
justified  because  at  some  point  the  decision  had  to  be  made 
whether  a  curve  was  acceptable  or  unacceptable. 
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TABLE  1.   Geometry  of  Locations  Studied. 
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RESULTS 


One  important  point  discovered  while  investigating  the 
problem  of  apparent  divergence  of  a  horizontal  curve  from  the 
tangent  was  that  when  the  observer  is  only  a  small  distance, 
such  as  3.5  feet,  above  the  plane  of  the  curve,  nothing  can  be 
done  to  make  the  roadway  appear  smooth.   The  tv;o  perspectives 
shown  in  Figure  6  illustrate  this  point.   Perspective  15  (Fig. 
6A)  is  a  one  degree  circular  curve  that  traverses  a  ten  degree 
deflection.   Perspective  31  (Fig.  6B)  is  the  same  location  but 
a  one-thousand  foot  spiral,  resulting  in  a  completely  spiralized 
curve,  was  used  in  place  of  the  circular  curve.   A  small 
difference  in  the  curves  can  be  detected,  but  neither  curve 
was  judged  to  be  visually  acceptable. 

However,  as  the  observer  is  raised  above  the  plane  of  the 
curve,  an  increasing  display  angle,  the  addition  of  spiral 
curves  to  the  horizontal  alignment  does  result  in  an  observable 
improvement  in  the  appearance  of  the  curve.   The  effect  of  dif- 
ferent spiral  curves  can  be  observed  in  Figures  7  and  8.   The 
sight  distance  and  display  angle  are  constant  for  all  perspec- 
tives.  The  length  of  spiral  curve  is  the  only  variable. 
Perspective  315  is  a  circular. curve  with  no  spiral.   Perspec- 
tive 323  has  a  250  foot  spiral  and  perspectives  315  and  331 
have  500  and  1000  foot  spirals,  respectively.   Note  the  increased 
smoothness  of  curve  from  perspectives  315  through  331. 

The  proposal  by  Pushkarev  (2)  that  a  completely  spiraled 
curve  will  appear  to  have  a  sharp  bend  in  the  middle  was 
investigated.   The  perspectives  in  Figures  7  and  8  illustrate 
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A.      No  Spiral    (Circular  Curve) 


B,      Completely  Spiralized 


FIGURE   6.        Effect  of   Spiral  Curves   at   Smll   Display 
Angles 
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FIGURE    7.         Effects   of    Spiral   Curves    (DA  ^- 
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B.      Spiral  Length   =   1000   Feet    (Coirpletely 
.    Spira lifted) 


FIGUPvT:    8.        Effects   of    Spiral   Curves    (DA  =    .015   Radieais, 
SD  =-■   1000   Feet,    D  =    1°    and  A   =   10°) 

Plate   2 
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that,  contrairy  to  appearing  sharp,  the  completely  spiraled 
curve  appears  smoothest  of  all  the  curves.   Therefore,  it 
may  be  concluded  that  the  length  of  spiral  curve  needed  at 
any  given  location  is  not  related  to  the  ratio  of  spiral  curve 
to  circular  arc.   There  are  tv;o  cases,  however,  which  were 
not  investigated.   The  first  is  when  a  vertical  curve  is 
superimposed  on  the  spiral,  a  problem  beyond  the  scope  of  this 
research,  and  the  case  where  the  observer  is  located  on  the 
curve  and  views  the  junction  of  the  two  spirals. 

An  investigation  was  also  undertaken  to  determine  if 
Godin's  recommendation  of  a  spiral  length  of  R/9  was  valid. 
As  was  pointed  out  in  the  previous  paragraph,  the  appearance 
of  any  given  curve  improved  with  additional  length  of  spiral 
added.   The  recommendation  to  use  a  length  of  spiral  equal 
to  R/9  is  further  limited  because  no  mention  is  made  of  the 
display  angle.   However,  for  almost  all  sight  distances  and 
display  angles,  it  was  felt  that  this  recommendation  resulted 
in  a  spiral  length  that  was  not  sufficient  to  be  visually 
significant. 

The  effects  of  various  display  angles  were  studied  to 
determine  the  extent  that  they  affect  the  appearance  of  a 
curve.   Figures  9  and  10  pictorially  show  how  the  display 
does  affect  the  appearance  of  a  curve.   All  perspectives  had 
the  same  geometry  and  sight  distance.   The  display  angles 
for  each  curve  are  listed  with  the  figure.   Note  that  as 
the  display  angle  increases,  the  appearance  of  the  curve 
improves. 
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h.      Display  Angle   -    ,007   Radian 


B.      Display  Angle   ~    .02  7   Radian 


FIGURE   9,         Effects    of   Display   Angle    (D  =^   2°,    A 

1      =   500    feet,    SD  =   500    feet) 
s 

Plate    1 


=   10°, 
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A,      Display "Angle   =    .06  7   Radian 


B,      Display   Angle   =    ,107    Radian 

FIGURE    10.      Effects    of   Display  ?j\qle    (D  ^-   2«,    A   =    10' 
1      =   500    feet,    SU   =  -SOO    feet) 
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An  attempt  was  made  to  construct  a  graph  showing  the 
relationship  bet\7een  sight  distance,  display  angle  and  the 
geometry  of  the  curve.   However,  the  wide  range  of  variables 
used  made  it  extremely  difficult  to  formalize  any  sort  of 
graph  including  all  of  the  varia±)les.   Figure  11  is  a  graph 
illustrating  the  relationship  between  sight  distance  and 
length  of  spiral  when  the  display  angle  and  degree  of  curve 
are  held  constant.   The  display  angle  and  degree  of  curve  for 
this  graph  were  chosen  so  that  the  number  of  available  data 
points  used  would  be  a  maximum,   A  quadratic  curve  was  chosen 
because  it  was  felt  that  at  long  sight  distances  no  spiral 
would  be  visually  acceptable. 
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RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

Late  in  the  research,  v/hile  discussing  the  problem  of 
determining  what  gives  any  given  curve  a  bad  appearance  with 
my  major  professor,  Dr.  Bob  L«  Smith,  it  was  brought  out  that 
the  amount  of  curvature  or  sharpness  of  curvature  of  a  curve 
in  the  perspective  picture  plane  should  be  a  good  indicator 
of  the  appearance  of  that  curve.   After  all,  that  is  exactly 
what  had  been  done  to  evaluate  the  perspectives,  i.e.  visually 
measuring  the  sharpness  of  the  curves.   It  was  decided,  there- 
fore, to  attempt  to  use  the  change  of  slope  in  the  perspective 
picture  plane  as  a  means  of  predicting  the  appearance  of  the 
curves.   Subsequent  discussion  made  it  seem  that  the  rate  of 
change  of  slope  of  the  curves  in  the  perspective  picture  plane 
would  be  of  more  value, 

A  Fortran  computer  program  was  prepared  that  would  convert 
the  points  of  a  curve  into  the  perspective  picture  plane 
coordinates  and  calculate  tl\e   slope  between  each  of  the  adjacent 
points.   The  change  of  slope  between  adjacent  slopes  was  then 
calculated  and  the  rate  of  change  of  slope  obtained  by  divid- 
ing the  change  of  slope  by  the  average  picture  plane  distance 
of  the  tv/o  slopes.   From  the  output,  change  of  slope  and  rate 
of  change  of  slope,  it  v/as  hoped  to  obtain  a  rate  of  change 
of  slope,  such  that,  any  rate  of  change  greater  than  the 
critical  rate  would  not  have  an  acceptable  appearance.   Only 
the  maxim.um  rate  of  change  for  each  curve  was  investigated 
as  this  was  thought  to  be  the  critical  point,  visually,  of  the 
curve. 
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The  initial  trial  of  this  hypothesis  resulted  in  unantici- 
pated results.   For  any  given  sight  distance  and  curve,  the 
maximum  rate-of-change  of  slope  decreased  as  the  appearance 
of  the  curves  improved,  as  was  expected.   However,  for  a  given 
viewing  angle,  the  relationship  did  not  give  the  anticipated 
result.   As  the  sight  distance  increased,  the  rate  of  change 
of  the  slope  decreased  while  the  appearance  of  the  curve  was 
becoming  poorer.   Inspection  of  the  change  of  slopes  revealed 
they  did  give  the  relationships  which  were  expected  from  the 
rate  of  change  of  slope. 

After  further  study,  the  situation  was  explained  by  the 
fact  that  the  points  were  all  located  at  equal  spacings  in 
the  space  {x,y,z)  coordinate  system  and,  therefore,  the 
changes  of  slope  were  actually  rates  of  change.   The  calculated 
rates  of  change  from  the  computer  program  were  calculated  from 
the  picture  plane  distance.   Although  the  points  were  equal 
distance  apart  in  the  space  coordinate  system,  the  coordinate 
transformation  into  a  perspective  view  made  the  picture  plane 
distance  of  the  points  near  the  observer  greater  than  the 
distance  of  the  points  further  from  the  observer.   There- 
fore, the  changes  of  slope  of , the  near  points  were  divided 
by  a  larger  number  than  the  more  distantly  located  points. 
However,  the  nature  of  a  perspective  picture  is  such  that  an 
observer  sees  all  the  points  as  being  separated  by  an  equal 
distance  when,  in  fact,  this  is  not  the  physical  case  in  the 
picture  plane. 
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Figure  12  is  an  illustration  showing  the  logic  of  using 
the  rate  of  change  of  slope  to  determine  the  visual  appearance 
of  a  curve.   Slope  lines  v/ere  drawn  through  corresponding  points 
of  each  perspective.   It  can  be  observed,  however,  that  the 
change  of  slope  is  greater  for  the  unacceptable  curve  than 
for  the  acceptable  one.   By  increasing  the  frequency  of  the 
calculation  of  the  slope  lines  to  the  extent  that  every  defined 
point  is  the  end  point  of  a  slope  line,  the  rate  of  change  of 
slope  can  be  calculated  for  each  point  and  the  maximum  value 
selected  for  comparison  to  the  critical  value. 

Table  2  shows  the  values  obtained  when  this  approach  v;as 
tested.   All  of  the  calculations  were  made  using  the  geometry 
of  one  curve;  D  =  30  minutes,  A  =  20  degrees  and  length  of 
spiral  =  2000  feet.   It  can  be  noted  that  the  numerical  values 
obtained  were  extremely  well  segregated  into  the  three  visual 
classifications  used  for  rating  these  curves. 

This  approach  is  not  restricted  to  points  which  are  equally 
spaced  in  three  dimensions,  hov.'ever.   The  complexity  of  the 
coordinate  transformation  requires  the  use  of  a  digital 
computer  and  it  is  a  small  matter  for  a  computer  to  calculate 
the  distance  betv/een  the  points  in  three  dimensions.   By 
dividing  the  change  of  slopes  by  half  the  total  distance  of 
the  two  slope  lines,  a  rate  of  change  of  slope  can  be  calculated. 
If  this  rate  of  change  is  less  than  the  critical  value,  the 
curve  will  appear  smooth  v/hen  constructed  in  three  dimensions. 
Conversely,  if  the  rate  of  change  is  greater  than  the  critical 
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FIGURE    32.      Illustration   of   Rate   of   Change   of   Slope 
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TABLE  2,   Rates  of  Change  of  Slope  Versus  Curve  Ratings 

Acceptable        Questionable        Not  Acceptable 

2.8  3.1  10.6 

1.6  3.4  11.7 

1.8  2.2  13.2 

2.0  3.8  14.6 

1.1 

1.2 

1.4 

1.5 

.7 
1.0 

.7 
1.0 

.7 
1.2 

.8 
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value  /   the   three-dimensional  curve  v;ill  appear  sharp   and 
"jerk"    avi/ay  from  the  tangent. 

The  preliminary   investigation  was  undertaken  with   the 
observer  stationed  three   feet  to  the   right  and  3.5    feet   above 
the   centerline  with   the    centerline  used  as  the   line   from  which 
the   slopes  were   calculated.      Use  of   any  other   line,   such   as 
one   of   the  edge  of  pavement   lines ,  would  be  valid  provided 
the   observer  is    located   in   the   same   relative   position.      The 
center  of   interest  was   on  the   centerline   approximately   at  the 
beginning   of  the   curve.      The   location   of   the   center  of   interest 
is  not  as   critical   to  the   analysis   of  rate  of   change   of   slope 
as   is   the   observer  position. 

In   recommending  this   particular   area  of   study   for 
additional   research,    it  should  be  pointed  out  that   the  basic 
concept  of   using  the   rate  of   change   of  slope   in  the  picture 
plane   seems   to  be  very   logical   and   straight-forward.      Whether 
the   formulation  presented  in   this   thesis   is  entirely  valid  or 
not  can   only  be   proven  by  additional  testing.      However,    it   is 
felt  that  this   method   for  determining  the  visual   acceptance 
of   a  curve  holds   great  promise. 
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CONCLUSION 

From  the  limitations  imposed  by  the  scope  of  this  study 
and  from  the  data  collected,  the  following  was  concluded: 

1.  Spiral  curves  do  improve  the  appearance  of  most 
circular  curves, 

2.  VThen  the  observer  is  near  the  plane  of  the  curve, 
there  is  no  significant  difference  in  the  appearance  of  a 
spiraled  and  an  unspiraled  curve, 

3.  As  the  distance  from  the  curve  to  the  observer  in- 
creases, the  length  of  spiral  needed  increases  for  good  visual 
quality, 

4.  As  the  height  of  the  observer  raises  above  the 
plane  of  the  curve,  increasing  display  angle,  the  length  of 
spiral  needed  decreases  for  good  visual  quality, 

5.  Curves  which  consist  entirely  of  spiral  curves  give 
the  best  visual  appearance,  all  other  conditions  being  equal, 

6.  The  rate  at  which  a  curve  visually  appears  to  diverge 
from  the  tangent  affects  the  visual  quality  of  the  curve. 
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IF(U I  5.5,10 
NSKIP(I)=NSKIP(  IX-1 

H(I,J)  =  0.   _  .  . 

V(I,J)=0. 

GO  TO  20 

H(I,J)=  SCALE*A2/(B*D)<-(XE*YY-YE*XX) 

V(I,J)  =  SCALE*A/ID*D)^MB2«ZZ-ZE*(XE*XX  +  YE*YY)I 

CONTINUE 

IF(NSKIP{18)-1)  60,60,50  

WRITE  (3,6^2) 

DO  '»0  1  =  1,  LI  N't 

NN=NSKIP( I )+l 

WRITE  (3,632)  { I , J,H( I , J ) , V( I , J ) , J=NN,N) 

PLOTN^NPLOT 

CALL  NUMBER  (+5 . , 0 . , 1 . , PL0TN,90. 0,-1 ) 

CALL  DRAW  ( L I NE , N , H, V ,NSK I P ) 
CALL  PLOT  (12. ,0., -3) 

NOPS=NOPS-I 

ir(NOPS)  30,30,3 
L0C=LCC-1 

IF(LOC)  898,898,1      ; 

CALL  PLOT  (  0.,0.,999) 
WRITE  (3,638) 

STOP  „^  

END 


40 
126  SUBRQUriNE  CLGK  ( XX, YY, ZZtNu, KEY, I SK IP ) 

MAIN  SUBROUTINE  FOR  CENTERLINE  COORDINATES 


C 

C 
C 

127 

900 

128 

901 

129 

902 

130 

903 

A31 

904 

C 
C 

132 

925 

133 

926 

134 

940 

135 

945 

136 

949 

137 

950 

138 

951 

139 

952 

140 

953 

141 

954 

142 

955 

143 

956 

144 

957 

145 

958 

146 

959 

147 

960 

148 

961 

149 

962 

150 

963 

151 

964 

152 

974 

153 

974 

154 

975 

155 

C 

976 

C 

156  " 

157 

158 

159 

C 
C 
C 

-  ■    - 

160 

161 

I 

162 

163 

110 

C 

* 

164 

c 

* 

165 

166 

167 

168.. 

READ  FORMA fS 
F0RMAT(10I5) 

FORMAT  (6F1 0.2)     _^ 

F0RMAT(I3,F10.2) 

rORMATC JP5.0,F3.0,F6.2) 

FORMAT{F5.0,F3.0,F6.2,6X,F5.0,Fi.0,F6.2) 


WRITE  FORMATS 

FORMAT! 53X,31HEM0  OF  ROADWAY  GEOMETRY  PROGRAM//) 

FORMAT  (34X,  70H*  *********«**********<*****************♦***********■< 

FORMAT  (IHl)       _       ^^^        __     

FORMAT!/////) 

FORMATdHl/////////) 

F0RMAT(10X,49H**«  ALIGNMENT  SIZE  EXCEEDS  DIMENSIONED  STORAGE  BY^i; 
1) 

F0RMAT(5X,54HWARNING  ♦**  NUMBER  OF  POINTS  MAY  EXCEED  DIMENSION  Sli 
IE) 

FORMAT(10X,62HTHE  TANGENT  DISTA^JCE  BETWEEN  TWO  SUCCESSIVE  CURVES  i 
IS  NEGATIVE) 

F0RMAT(1M+,100X,17HTANGENT  D  IRECT ION, F14.2 ) 

FORMATdHl, 27X,25HHORIZONTAL  CURVE  GEOMETRY/////) 

F0RMAT(28X,23HVERTICAL  CURVE  GEOMETRY/////) 

FORMATdHl, 47X,9HTHERE  ARE,  I3,24H  POINTS  IN  THE  ALIGNMENT) 

FORMAT(50X,22I1STATION  OF  FIRST  POINT,F10.2) 

FORMAT150X,21HSTATION  OF  LAST  P0INT,F11.2) 

FORMATd0X,5dH*),38H  INVALID  CROSS  SECTION  INDICATOR  VALUE  ) 

FORMAT (36X,5HCURVE, 13) 

FORMAT! 36 X,8H*** ****«) 

FORMATdHl, 30X,  19HCENTERLINE  GEOMETRY/) 

F0RMAT{17X,I5,4F10.2) 
F0RMAT(21X,1HI,6X,1HX,9X,1HY,9X,1HZ,8X,3HSTA) 

FORMAT!//////////)  _         _ 

F0RMATd01X,23HDISrANCE  BETWEEN  CURVES, F8. 2) 
FORMAT !101X,14H TANGENT  LENGTH, Fl 7, 2 ) 

F0RMAT!//48X,23HCR0SS  LINES  DRAWN  EVERY, F8.0, 5H  FEET/48X, 24HP0INT: 
1  ARE  DEFINED  EVERY, F7.0, 5H  FEET) 

COMMON  AND  DIMENSION  STATEMENTS  

DIMENSION  XX!15, 90), YY!15, 90), ZZ!15, 90), ZdOO)    ' 
DIMENSION  SPIRdO),DEFL!10),TANL!ll) 

COMM0N/CL/XI90) ,Y!90) ,5TA!100)  

COMMON/RDl/  OIST,I,N,NDEFL!10),RDEFLd0),0Cd0)^ 


READ  DATA  AND  INITIALIZE  VALUES 
DO  1  1=1,10 

SPIR!I)  =  0.  _       _   __.    

DO  110  1=1,60 
STA!I)=0.0 

READ  NUMBER  OF  PIS!NPHI)  AND  PVISINPVI)    _    __   „ 

READ  11,900)  NPHI,NSPIR,NPVI,KEY 

READ  DISTANCE  BETWEEN  POINTS  AND  INITIAL  DIRECTION  FROM  NORTH 

READ  !1,903)  DI  ST,  XDI  ST,  A,  B  ,C  __  _     

THETA=!A  +  B/60.  +  C/3600.  )*0. 0174533  '"'     '         "' 

NTAN=NPHI+1 

IFIXDIST-DIST)  109,117,115 


169 

109 

170 

117 

171 

172 

115 

173 

174 

116 

175 

100 

176 

177 

178 

101 

179 

106 

180 

181 

C 

* 

182 

105 

C 

c 

c 

♦ 

183 

184 

185 

186 

112 

187 

188 

113 

189 

114 

190 

107 

c 

♦ 

191 

192 

193 

194 

102 

195 

196 

197 

103 

198 

111 

199  ^ 

200 

201 

202 

108 

i03 

204 

104 

^205 

c 

* 

205 

c 

211 

c 

c 

„  207 

210 

208 

209 

210 

211 

212 

213 

214 

99 

215 

216 

217 

c 

97 



._.  c 

41 
XDIST=DIST 

ISKIP=1 

GO  TO  116         ._  

ISKIP=XDIST/DISr 
XDIST=DIST*FLOAT(ISKIP) 

IF(10-NPHI)100,10l,101   

IXCES=NPHI-10 
WRITE  (3,950)  IXCES 

GO  TO  899 

IF(8-NPVI)  106,105,105 
IXCES=NPVI-8 

WRITE  (3,9501  IXCES   

GO  TO  899 

READ  COORDINATES  OF  INITIAL  POINT  AND  STATION 

READ  (1,901)  Xd),  Yd),  ST  Ad)  

DIRECTION  CHANGES  ARE  READ  AS  DEFLECTIONS 

CLOCKWISE  IS  POSITIVE  AND  COUNTERCLOCKWISE  IS  NEGATIVE 

READ  DEFLECTION  ANGLES  (DEFL)  AND  DEGREES  OF  CURVE  (DC) 

DO  107  I=1,NPHI 

READ  (1,904)  A,B,C,D,E,F 

IF(A)  112,113,113      

DEFL(I)=A-B/60.-C/3600. 
GO  TO  114 

DEFL(  I)  =  A  +  B/60.+C/3600.  ^ 

DC{I)-D+E/60.+F/3600. 
CONTINUE 

READ  TANGENT  LENGTHS  . 

READ  (1,901)  (TANLd ),I=1,NTAN) 
TOTAL=0. 

DO  102  I=1,NTAN    

TaTAL=TOTAL*TANL{I ) 
NTOT=TOTAL/DIST 

IF(NT0T-90)lll,lll,103   

WRITE  (3,951) 
DO  104  I=1,NPHI 

RDEFL(I)=DEFL(I )*. 01745  33  

NDEFL{I)=1 

IF(DEFLd)  )108,104,104 

N0EFL(I)=-1  , 

DEFL(I)=-DEFLII ) 
CONTINUE 

IF(NSPIR)2Il,2lO,2il  .; 

READ  SPIRAL  CURVE  NUMBER  AND  LENGTH 
READ  (1,902)  (N,SPIR{N),I=1,NSPIR) 
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TD1=0. 

N=2 

S1=0.        ^  _ 

1  =  1 

WRITE  (3,954) 

GO  TO  97  _  _ 
99  1=1+1 

DHETA-THETA*57.3 

WRITE  (3,953)  DHETA 
97  IF(SPIRd)  )3,3,4 


PRELIMINARY  CIRCULAR  CURVE  CALCULATIONS 


218 

3 

219 

20 

220 

221 

22 

222 

23 

223 

224 

225 

226 

227 

209 

228 

229 

2 

230 

27 

231 

232 

233 

234 

6 

235 

8 

236 

237 

5 

238 

239 

240 

241 

242 

17 

243 

244 

245 

7 

246 

247 

9 

248 

249 

250 

251 

252 

18 

253 

254 

255 

256 

257 

258 

259 

c 

c 

260 

4 

261 

24 

262 

263 

25 

C 

264 

265 

266 

267 

268 

269 

270 

26 

271 

272 

273 

274 

42 
IF{NTAN-I)20,20,22 
TD2=0. 

GO  ro  23  _  _._ 

TU2= 5729.58* TAN {.0087267*DEFL( I) )/UC(I) 
TT=T2 

TB  =  TANL(I  )-TDl-TD2 _ 

WRITE  (3,974)  TB 
WRITE  (3,975)  TD2 

IF(TB+10. )209,2,2   .   

WRITE  (3,952) 
GO  TO  899 

IF(T2)6,6,27  _ 

X{N)=PTX+T2*SIN(THETA) 
Y(N)=PTY+T2*C0S(THETA) 

STA(N)  =  PTS  +  T2  _   _ 

N=N  +  1 
IF(TB-DIST)17,8t8 

TT=TT+DIST  _       

IF(TT-TB)5,7,7 
X(N)=X(N-1)+DIST*SIN(THETA) 

Y(N)=Y(N-1)+0IST*C0S(THETA)   

STA(N)  =  STA(N-1)  <-DIST 
N=N+1 

GO  TO  8        _   __ 

Cl=OIST-TB 
Dl  =  TB 
GO  TO  9 
C1=TT-TB 
D1=0IST-C1 

ANGLE=THETA  _  .  

PCX=X(N-1)+D1*SIN(THETA) 
PCY=Y(N-1)+D1*C0S(THETA) 

PCS=STA(N-1)*^D1 

IF(I-NPHI)18,18,799 
WRITE  (3,940) 

WRITE  (3,945)     _  

WRITE  (3,960)  I 

WRITE  (3,961) 

CALL  CIRCLE  (CI , ANGLE, T2 , DEFL ( I ) , PCX, PCY, PTX, PTY, PCS, PTS ) 

THETA=THETA+RDEFL(I) 

TD1=TD2 

GO  TO  99 


PRELIMINARY  SPIRAL  CURVE  CALCULATIONS 
IF(NTAN-I )24,24,25 
T02=0. 
GO  TO  26 

RC=5729.58/DC(I  )  

THES  IS  IN  RADIANS 

THES=SPIR( I)/(2.*RC) 

YS=SPIR( I )*THES/3, 

XS=SPIR{I )-(SPIR( I )*THES**2/10.  ) 

P=YS-RC*(1.-C0S(THES) ) 

OK=XS-RC*SIN(THES) 

T02=(RC+P)*TAN( . 008  7267*0EFL ( I ) )+DK 

TB=TANL(I)-TD1-T02 

WRITE  (3,974)  TB  _ 

WRITE  (3,975)  TD2 

TT=T2 

IF(TB+10. )10,11,11 


292  _ 

15 

293 

294 

14 

„295  „ 

16 

296 

297 

19 

298  „_ 
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275  10  WRITE  (3,952) 

276  GO  TO  899 

277  11  IFIT2)2R,28,29         

27'3       29  XCN)  =  PTX+T2*Sirv!{TH£rA) 

279  Y(N)=PTY+T2*C0S(THETA) 

280 STA(N)  =  PTS  +  T2  

281  N=N+1 

282  28  IF(TB-DISTJ21, 12,12 

283 12  TT=TT  +  DIST  . . 

284  IF(TT-TBH3,14,15 

285  13  XlN)=X(N-l)+DIST*SIN(TnerA) 

286  Y{N)=Y(N-1)+DIST*C0S( THETA)   ; 

287  STA(N)=STA(N-l)+Dr5T 

283  N=N+1 

28'J  GO  TO  12 : 

290  21  S1=DIST-TB 

291  GO  TO  16 

GO  TO  16 
S1  =  0. 

AMGLE=THETA        1_ 

IF(I-NPHI)19,19,799 

WRITE  (3,940) 

WRITE  (3,960)  I  

299  WRITE  (3,961) 

300  CALL  SPIRAL  ( SI, ANGLE, T2, SPIR, rHES,PTX,PTY,PTS) 

301  THETA=THETA4-RDEFL(  I  ) 

302  TDl=TD2 

303  GO  TO  99 

304  799  N=N-1        .„  _  

305  WRITE  (3,956)  N 

306  WRITE  (3,957)  STA(l) 

307   WRITE  (3,958)  STA(N)   

308  WRITE  (3,976)  XOIST,DIST 

309  WRITE  (3,940) 

310  WRITE  (3,955)  .  _  .  

C    *  SUBROUTINE  VERT  CONTAINS  READ  STATEMENTS 

311  CALL  VERT  ( NPVI , N, STA, 7 ) 

312  ._ WRITE  (3,962)  

313  WRITE  (3,964) 

314  DO  30  1=1, N 
315 30  WRITE  (3,963)  I  ,  X  (  I ) , Y( I ) , Z ( I ) , STA(  I ) 

316  IF  (KEY)  31,32,33 

317  31  WRITE  (3,959) 

318  _  _    WRITE  (2,963)  N  

319  WRITE  (2,963)  (  I , X ( I ) , Y ( I ) , Z ( I ) , STA ( I ) , I^l,N) 

320  GO  TO  899 
321 32  J  =  l        _. 

322  GO  TO  34 

323  33  IFCKEY-l)  32,32,35 

324  35  J  =  9  _  „   

325  34  DO  920  1  =  1, M 

326  XX{J,I)=X(I) 

327  YY(J,I)  =  Y{I)         ^   ^    ,   „_    _   ^ 

328  920  ZZ(J,  I)  =  7.(  I) 

329  999  WRITE  (3,949) 

330  ^ WRITE  (3,926)     

331  WRITE  (3,925) 

332  WRITE  (3,926) 
___333  L WRITE  (3,940)    _ _.  .  _ 
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33't 

NrNj=N 

3  3? 

RETURN 

336 

899 

MN=0 

337 

RETURN 

338 

END 

.  339 

1 

3A0 

1 

:_3^i    _„ 

1 

342 

.   343 

10  ' 

344 

1 

>  345 

! 

346 

1 

347 

1 

348 

.   349 

_  350 

^  351 

1 

.  352 

1 

.  .353.    _ 

354 

( 

.   355 

356 

40  ' 

-  357 

1 

.  358 

359 

360 

1 

,      361 

362 

„_50 

^  363 

,  364 

1 

_365 

•   366 

.   367 

60 

.  368 

369 

SUoRObTINE  0FRD2  ( X, Y t Z, NPPL J 
DIMENSION  X(15,90),Y(15,90),Z(13,90) 

DO  10  I=I,NPPL     .       

X(9,I )=X{l,I) 
Y(9, I  )  =  Y{1,I) 

NLINE=5         _.  .  _._ 

NM1=NPPL-1 

DO  50  K=2tNLlN£ 

D=-22.0  

IF(K.EC.3)  D=-12.0 
IF(K.EQ.4)  0=12.0 

IF(K,EQ.5)  0=22.0  

DO  40  1=2, NMl 
L=I+1 

J=I-1  

CALL  DELTA  ( X,Y , L, J ,D,CX,CY ) 
X{Kt I )=X{l,I )+CX 

Y(K,I )=Y(l,I )+CY  

CALL  DELTA  ( X , Y , 2 , 1 , D ,CX , CY ) 
X(K,1)=X(1,1)+CX 

Y(K,1)=Y{1,1)+CY  

CALL  DELTA  ( X, Y, NPPL, NMl, D,CX ,CY ) 
X(K,NPPL)=X{ 1,NPPL)+CX 

Y(K,NPPL)=Y(1,NPPL)+CY  _ 

DO  60  K=2.N'LINE 
CZ=1.0 

IF(K.EC.3.nR.K.EQ.4)  CZ=0.25 

DO  60  1=1, NPPL 

z(K,  I  )=z(i,n-cz 

RETURN  

END 
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370 

371 
372 
373 
374 
375 
376 
377 
378 
379 
380 


SUBROUTI^IE  DELTA  (  A  ,B  ,M,L  ,01  ♦CX,  CY) 

CALCULATE  COORDINATE  CHANGES  FOR  ALL  POINTS  FROM  CENTERLING 

DIMENSION  A(15,90) ,B(15,90)  

IF(A{9,M)-A(9,L) )30,3l,30 

31  THFTA=1.5708 

GO  TO  32  

30  ANGLE=ABS( ( B ( 9, M ) -B ( 9,L ) ) / I A( 9, M )-A( 9, L) ) ) 

THETA=ATAN{ANGLE) 

32  CX=DI*SIN(THETA)         __. 

CY=DI«COS(THETA) 

RETURN 

END  


47 

381  SUBROUTINE  0FS-N2  (X,Y,2,N'PPL) 

C      OFFSET  SUBROUTINE  FOR  A  TWO  LANE  ROADWAY 

DIMENSION  ST  A  ( 90  )  ,  X  ( I  ^ ,  90  )  ,  Y  ( 1 5 ,  90  )  ,  Z  ( 15 ,  90  ) ,  DR(90 )_,  DL  190] 

DIMENSION  BSR(90) ,BSL{90) 
F0RMAT(4F10.3) 

NM1=NPPL-1  .   

CUT  {-)      FILL  1+)    FROM  CENTERLINE  OF  ROADWAY 

DEPTH  OF  CUT  OR  FILL  ON  RI3HTIDR)  AND  LEFT(DL)  SIDE  OF  ROADWAY 

BACKSLOPE  ON  RIGHT(BSR)  AND  LEFTIBSL)  SIDE  OF  ROADWAY 

READ  DEPTHS  AND  BACKSLOPES 

READ  (1,102)  (DR(I),BSR(I) ,DL( I ) ,BSL(I), I=1,NPPL) 

WRITE  (3,102)  DR  (  NPPL  )  ,  BSR  (  NPPL  )  ,  DL  (NPPL  >  ,  BSL  (  MPPL  )_ 

CALCULATIONS 
NLINE=6 

DO  50  K=1,NLINE 

IF(K.EQ.2)  D— 22. 
IF(K.EG.3)  D=-12. 

I  F  {  K .  EQ .  't )  D=  12  . 

IF(K.EC.5)  D=22. 
DO  AO  1=2, NMl 

IF(K.EQ.l)  D=-22.-ABS{DL{I)*BSL( I ) )   

IF{K.EQ.6)  D=22.+ABS(DR(I )*DSR(I ) ) 
L=I  +  1 

J=I-1  

CALL  DELTA  ( X,Y, L , J , D,CX,CY ) 
X(K,I)=X(9,I )+CX 

Y(K, I )=Y{9, I )+CY  

CALL  DELTA  ( X, Y , 2 , 1 , D ,CX, CY ) 
X(K,1)=X(9,1)+CX 

Y{K,1)=Y(9,1)+CY  __  

CALL  DELTA  { X ,Y , NPPL , NMl , D,CX, CY) 
X(K,NPPL)=X{9,NPPL)+CX 

Y(K,NPPL)=Y(9,NPPL)+CY   

DO  60  K=1,NLINE 
00  60  1=1, NPPL 

IF(K.EQ.l)  CZ=DL(I)         _  

IF(K.EQ.2  .0R.K.EQ.5)  CZ=1.0 
IFIK.EQ.3  .OR.K.EQ.A)  CZ=.25 

IF(K.EQ.6)  CZ=DR(I)    

Z(K,I )=Z(9,I)-CZ 

RETURN 

END : 
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384 
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C 

C 

C 

C 

« 
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C 
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60 
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!.  "Vri,    i ,    rr;;-)   3 
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DATE    =    69162 
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16/^o/^ 


SUBKOUTINt.    ni'LN     (  X  ,  Y  ,  7  ,  S  T  A  ,  NP  ,Kb"Y  ) 

CP'CriSIO  J    X(  ir>,300),Y(  15,5:0)  ,Z{  l-i.iJOI.SlAt  i:rj),  STAT  (2,  IOC) 

C.I  MP  ?.'!">  I  Ci    [-1>V(?,  130)  ,THr.TA(  jj.j) 

i\n^M"-i ■ 

If  (K^Y.r.T.A  )     'JC;    T'.i    101 
GO    TP.     (1,  1.11,  Id,;?)  ,KEY 

1  cn  =  27. 

NCL^  1 
CD    1"     ^ 

2  CL=?'t. 

3  cfi   3.    1=1 ,  ;n 

Ih(XC:Cl  ,  I  M  )-X(^;CL,  I  )  )    2:,1;',20 
iC    Tl,!  I  .'  (  I  )  ^  1.5  7.B 

GO    TC     3r 
?u    ftV'M^--  '.'•.,  (  (Y(MCL,  1*1  )-Y(.\LL,  I  )  )/(X(-JLL,It  1  )-/.  ("ICL,  I)  )  ) 

Th:  I  ;.  (  I  )  =  AI  «  ,(  A.\3LEi 

?'j  cr  \'i  I  -'ur 

cr   '.''    !=/,Nn 
/.u  Tiff  t;.  (  1 1^ :  TH:  TA(  i-i)MHi:rA(  I )  1/2. 

C  *    flLAf;    -iOM^LK    ilF     I^!FL^CWCJ\     PHlXTi,    Fi'il    r^I'.HT     A  jO    I  i  K  T    i;  CK  SL  Ci'c' S 

Ri-AC    (i,95C)    \':-<,f:L 
c       *  Ri  A[;  SJArio"4  A'^o  Lirv  vT  u: .  fo^  ^>riiir  a-;u  Lf-t-i   qacksli;!'^: s 
f<L"*i'   (i,>>n    (siAi  ( 1 , 1)  ,n  =  v(i ,  I ) ,  1  =  1  ,m;<) 

K['P     (l,'^:')!)     (STAI  (2,  n  ,ri.rV(  2,  I  )  ,  1  =  1,     L) 

ES^6. 

DC    1::    K=l,2 

SI^»M^»  1.  J 

N=0 

U  (KL-Y.C'J.A)     •j=12 

NCL-\U 

I'   (K'l)     S),'..b,'i' 

\=1  . 

sin--'- --1.0 
5'>   uLi    !;■:    1  =  1.,  .^' 

LO    6'.      1  =  2  ,^  -i 

[F(S1  A(  n-  ■>!  -.  [  1K,..I)  '      7'  ,  c3:),6- 

60  cr-;T  r;^;*: 

Kt  TU;''i 

70    7CI,  n  ■ -■  L  "7(:<,  :     i  ;  »■(  SIA(  I  )-  bTATCK,  J-1)  )i  (rL^.vMK,  J)-'^Lt-V<K,  J-n  )/(S 
ir.M  (,■.,.!)-•,!  M  !-,..-l  )  ) 

cisT^.v +,-;-; -^^..^f  '•     .,  n-MNCLi  I  )  ) 

Cf   jr   'jc 

PO     /('!,  I  )-^L  :V  (:.,.'  ) 

cisr  =  Dr; 

9C     Xr.,  l)=X(  nCL  ,  ;  :  •!".ir;v<v-)iS7^bI  .(  IHETAl  I  )  ) 

Y(^M  )->(  iCu ,  •  i4biONtDisi<^Cijo{  th:ta(  n  ) 

101     .••  ',  IT^     M,   7  ■;■  )     rvf  Y 

-j'V  rr!N'.r(->h  -      ^,  n,  iox,.3'^<Hr;v\Lio  valu-^  or  key   i  j  hfli  suL'^Muri  -c) 
i)-)  "-u.  ru.-  .' 

•^00    r:.M'  AT  CI  ■ ) 

■>01     r-^:  ••  AT  (^^  ■  "  ,.   )   ■ 


V 


N    IV   3   LEVEL    1,    HOD    3                               OFLN                                    DATE   =    69U2                        l6/Ai,M^ 
9A_0    FORMA  r(10X,«  ERROR    STATEMEMT    NO.    60     .OFLNJSUBROUTJNEl) 

■   ""end 


^1  IV  C  LEVr-L  1,  MOD  3 


^, 


OrOUR  DATt  =  6V162  l6/45A«j 


SUftROUTINE  OHCUR  (  A  ,B  ,  M,  L,  01  ,CX  ,CY  )    .     ,  ,^^.-v..,- 

CALCULATH  COORDIMATE  CHANGES  FOR  ALL  POINTS  FROM  CENTERLIMfc 

D I  MENS  I  CNJ  A  ( 1 5 ,  90  )  ,  B  { 1  5 ,  90  ) 

ir(A(13,N)-All;J,L))30,31,30  

31  THErA=1.^703 

CO  TO  32 
30  A^JGLF  =  AHS(  (B(13,r.)-R(13,L)  )/(A(l3,M)-A{13,L))) 

THETA=ATAN(A\GLF) 

32  CX  =  [)I'^SrMTHETA) 

CY  =  t)I*CijS(THZrA)  .  - 

RETURN 

FNO 


1 


H 


IV 


LF.VEL 


10 


U= SIGN* 54, 

n=SIGN*42. 

D=SIGN*3j. 
0=SI3N*2A. 


1,    MOD    3  OFAL 

SUf>RCUTTME    O^AL    (  X,  Y./'.NPPL  )  _      

DI  MENS  I  Of)    X  (I'j,  90),  Y{  15,90  ),Z(15«90» 

ASSUMES    A8    FT    MEDIAN,    6    FT    INSIDE    AND 

DO    10    I=1,NPPL  __ 

X{13,I  )  =  X(9, I) 

Y(13,n=Y(9,I) 

Z(13,I)-Z{9,! ) 

NLrJE=ll 

SI0N=-1. 

NMl=NPPL-l 

DCI    SO    K  =  ?,NLlMt 

IF(K.0F.7)     SICM=+1. 

C=SI&N*6A. 

1F(K,'  O.J.Oi^.K.EC.lO) 

IF(K,'  J.A.n,<.K,EQ.9) 

iFiK.^c.s.n.i.K.Eu.a) 

IF(K.rf..G.iJ'<.K.E0.7) 

DO    AC-     I=2,Nf-;l 

L=I»1  J-  

J=I-l 

CALL    nrCUii     (X,Y,L,  J,n,CX,CY) 

X(K, I )=X{ 13, I )+CX 
AG    Y(K,  n=Y(13,  I)+CY 

CALL    I1FCUK     (X,Y,2,  1,0,CX,CY) 

X(K,l)=X( 1^, l)+CX 

YtK, 1)=Y( I3,1)*CY 

CALL    DFCu'^     (X,Y,NPPL,NKl,l),CX,CY) 

X(K,\'Pl'H=Xl  13,-^JPPL)+CX  _ 

50    YIK,NPPL)=Y113,MPPL)*CY 

DCI    6C    K  =  2,MLINE 

CZ=1.44 
ir(K.F0.3.nR.K.Et).5.CR.K.EG.8.n;<.K.EC.I0) 

IFCK.EO.A.OR.K.EC.O)    CZ=J. 

IF (K. E 0.6. OR. K. EC. 7)     CZ=.9A  . 

CO    60    I  =  1,\'PPL. 

60    Z(K,  I  )=Z( 13,  I  )-CZ 

RETUH-I  _ -      ..  - 

END 


DATE    =    69162 


1&/A5^ 


1 


10  FT  OUTSIDE  SHOULDERS 


CZ  =  .  19 


r. 


^17 

410 

A19 

.  420 

951 

421 

952 

422 

953 

423 

954 

424 

955 

425^ 

956 

H  426 

95/ 

427 

950 

428 

959 

429 

430 

431 

-  432 

433 

8 

434 

„ 1 

435 

436 

5 

437 

4 

.  438- 

439 

440 

441 

442 

443  . 

3 

'  444 

2 

445 

6 

446  . 

447 

448 

9 

449 

.  _.  10 

450 

451 

'452 

C 

453 

454 

455 

.  456 

457 

458 

23 

459 

460  . 

24 

461 

25 

,  462 

463 

464 

465 

SUDROUTINE  CIRCLE  (CI , ANGLE , T2, DEFL, PCX, PCY,PTX, PTYv PCS, PT 
COHMON/RDIV  OIS T , I , N, NDcFL 1 10 ) » ROEFL ( 10) , DC ( 10 ) 

COMM0N/CL/X(90) ,Yi90),STA(100)       

F0RMAT(////31X, 14HCIRCULAR  CURVE/ ) 
F0RMAT(3IX,1HX,16X,1HY,15X,3HSTA/) 

FORMAT ( 16X , 2MPC , 3F 1 7. 2 )  

F0RMAT(16X,2HPT,3F17.2) 

FORMAT(//24X,23HDEFLECriON  TO  THE  RIGHT, F7.2,8H  DEGREES) 

FORMAT(//24X,22HDEFLECriOfM  TO  THE  LEFT,F7.2t8H  DEGREES) 

FORMAT(/32X,7HRADIUS=»F10.2) 
FOPsnAT{/30X,l7HDEGREE  OF  CURVE  =,F5.2) 

F0RHAT{/32X,7HLENGTH=,F10.2)  ,  

RAD-5729,58/DC{  I ) 
CUR=100.*DEFL/DCII ) 

TC=C1  _    _    

IF(CUR-IC)  6,8,8 
IF(Cl)  3,3,1 

DA=TC*DC( I )*. 000087267  

IF(NDEFL( I) )5,5,4 
DA=-DA 

BNGLE^ANGLE  +  DA  _ [ 

X(N)^PCX+TC*SIN(BNGLE) 
Y(N)=PCY+TC*CCS(BNGLE) 

STA(N)  =  PCS*TC    _ 

N=N+1 
IF(CUR-TC)6,3,3 

TC=TC  +  DIST       _ _  _. 

IF(CUR-TC)6,1,1 
T2=TC-CUR 

DA=CUR*DC( I)*. 000087267      . 

IF(NDEFL( I ) )9,9,10 
DA=-DA 

ANGLE  =  ANCLE  +  DA  .     _       _.    _   _ 

PrX=PCX+ SI N( ANGLE )*CUR 
PTY=PCY+CnS( ANGLE )*CUR 

PTS  =  PCS  +  CUP  .  ._  _.. 

WRITE  CURVE  GEOMETRY 
WRITE  (3,951) 

WRITE  (3,952)  __  ^_ _ 

WRITE  (3,953)  PCX,PCY,PCS 

WRITE  (3,954)  PTX,PTY,PTS 

IF(NDEFL(I) )23,23,24  _^_„^    ___ 

WRITE  (3,956)  DEFL 

GO  TO  25 

WRITE  (3,955)  DEFL 

WRITE  (3,958)  DC  (  I  ) 
WRITE  (3,957)  RAD 

WRITE  (3,959)  CUR  _ 

RETURN 
END 


... 

':65 

467 

468 

.  469 

470 

952 

_  471 

.  953 

472 

954 

473 

955 

474 

956 

,  475 

964 

476 

965 

477 

966 

478 

967 

479 

974 

480 

-  481 

482 

483_  . 

..   26 

C 

484 

27 

_  485  . 

1 

•  486 

487 

-  .488 

489 

490 

„491 

-  492 

493 

494 

2 

495 

496 

497 

498 

499 

.  500  _ 

5 

501 

4 

502 

503  ._ 

504 

,  505 

_  506 

507 

7 

508 

3 

509 

510 

8 

.  511 

512 

11 

513 

12 

514 

515 

516 

.  517 

518  _ 

519 

6 

520 

521 

522 

9 

r  523 

10 

_  524 

SUBROUTINE  SPIRAL  (  SI ,  Ar^GLE  ,  T2,  5P  IR,TMi:S»  5TX,STY,  STS) 
COMMON/RDl/  DIST, I , N, NDEFL (10) , RDEFL ( 10 ) t DC ( 10 ) 

COiMM0N/CL/X(90)  ,Y(90),STAC100)  

OIMENSIO^J  SPIR(  12) 

r-QRMAT(3lX,  IHX,  16X,  IHY,  15X,  3HSTA/  ) 

F0RMAT(16X,2HCS,3Fi7.2)  

FORMAT (16X,2HST,3F17.  2) 

F0ilMAT(////28X,  20HSPIRAL  CURVE  (AHEAD)/) 

F0RMAT(////28X,21HSPIRAL  CURVE  (BEHIND)/) 

F0RMAT(16X,2HTS,3F17.2) 
F0RMAT(//27X,25HLIMITING  DEGREE  OF  CURVE=,F5,2) 

F0RMAT(/30X,1AHSPIRAL  LENGTH=, F 10. 2  )       

FORMAT (/31X,9HTHETA(S)=,F5.2,8H  DEGREES) 
F0RMAT{16X,2HSC,3F17.2) 

DHES=THES*57.30 

DISTS=100. 
IF(DIST-100. )  26,27,27 

DISTS=DIST  

CALCULATIONS  FOR  SPIRAL  TO  CIRCULAR  CURVE 
IF(S1)1,1,2 

TSX=X(N-1)4^DIST*SIN(AMGLE)    

TSY  =  Y{N-l)+DISTi=COS(ANGLE) 
TSS=STA(N-1)+DIST 

X(N)  =  TSX   ...  

YIN)=TSY 
STA(N)=TSS 

N=N+1 ,       

TS=0. 

GO  TO  7 

T1=0IST-S1 

TSX=X(N-1)+T1*SIN(ANGLE) 

TSY=Y{N-l)+Tl=f=COS(  ANGLE) 

TSS=STA(N-1)+Tl  .     __ 

DA=(Sl/SPIR( I ))**2*THES/3. 

IF(NDEFL{ I ) )5,5,4 

DA=-DA 

ANGS=ANGLE+DA 

X(N)=TSX+S1*SIN(ANGS) 

Y(N)=TSY+S1*C0S{ANGS)    ^   

STA(N)=TSS+S1 
N=N+1 

TS=S1  

IF(SPIR(I)-TS)6,6,3 

TS=TS+DISTS 

IF(SPIR( I)-TS)6,6,8  _ 

DA=(TS/SPIR( I ))**2*THES/3. 

IF(ND£FL(I ) )11,11,12 

DA=-DA 

ANGS=ANGLE+DA 

X(N)  =  TSX4^TS*SIN(ANGS) 

Y(N)=TSY  +  TS*CCS(ANGS)    ._  .. 

STA{N)=T5S*TS 
N=N+1 

GO  TO  3  „     

C1=TS-SPIR(I  ) 

DA=THES/3. 

IF(NDEFL(I ) )9,9, 10  _ 

DA=-DA 

ANGS=ANGLE+OA 

SCX=TSX4-SPIR(I)*SIN(ANGS) 


54 

525  SCY-TSY  +  SPIR(n*CuS(ANGS) 

526  SCS=TSS+SPIR(n 

_ C_     WRITE  SPIRAL  AHEAD  GEOMETRY     . ^ 

WRITE  (3,955) 
WRITE  (3,952) 

WRITE  (3,96^)  TSX,TSY,TSS  

WRITE  (3,97A)  SCX,SCY,SCS 
WRITE(  3,965)  DC(I) 

WRITE  (3,966)  SPIR(I)      _     _  

WRITE  (3,967)  DHES 
IF(NDEFL(I ) )2l, 21,24 

ANGLE  =  ANGLE-THES   __  _  .   

GO  TO  25 

ANGLE=ANGLE+THES 

CDEL  IS  CENTRAL  CURVE  DEFLECTION  IN  DEGREES  (ALWAYS  POSITIVE} 

CDEL={ABS(RDEFL(I) ) -2 .*THES ) *57. 30 

ANGLA=ANGLE 

CALL  CIRCLE  (CI , ANGLA, S1,CDEL, SCX, SCY,CSX,CSY, SC S,CSS )  

COEL=CDEL«. 0174533 
IF{NDEFL(I) )19, 19,13 

CDEL=-CDEL  ■ ^ 

ANGLE=ANGLE+CDEL 

CALCULATIONS  FOR  CIRCULAR  CURVE  TO  TANGENT 

TS=S1  _       _     _        _. „ 

IF(S1)14,1A,15 

TS=TS+DISTS 

IF(TS-SPIR( I))13,20,20  __ 

OA=TS*DC{ I )*.C00087267-(TS/SPIR(I ) )**2*THES/3. 

IF(NDEFL(I ) )16, 16,17 

DA=-DA  

ANGS=ANGLE+DA 
X(N)=CSX+TS*SIN(ANGS) 

Y(N)=CSY+TS*COS(ANGS)     _ 

STA(N)=CSS+TS 

N=N  +  1 

GO  TO  14  _  ^   _  _  _  _     _   __  ^    __   _ 

T2=TS-SPIR(I  ) 
DA=2.*THES/3. 

IF(NDEFL(I  ))  22,22,23 ^ 

DA=-DA 
ANGS=ANGLE+DA 

STX=CSX+SPIR(I)*SIN(ANGS)       _  _   __ 

STY=CSY+SPIR(I)*COS(ANGS) 
STS=CSS+SPIR(I) 

WRITE  SPIRAL  BEHIND  GEOMETRY  __  

WRITE  (3,956) 
WRITE  (3,952) 

WRITE  (3,953)  CSX,CSY,CSS  __ 

WRITE  (3,954)  STX,STY,STS 
WRITE  (3,965)  DC (  I) 

WRITE  (3,966)  SPIR(I) 

WRITE  (3,967)  DHES 

RETURN 

END 
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►  533 
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21 
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24 

C 
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25 
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19 

1  544 

13 

C 
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547 

14 
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15 
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18 

550 

„  551  „ 

16 

552 

17 
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554 
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.  556 

557 
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20 

559 

560 
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22 

>  562 

23 

563 

564 

565 

C 

566 

f  567 

568 

569 

570  . 

571 

572 

573 

„  574 

575 

SUBROUTINE  VERT  1 NPVI  .NNN.STA.ELEV)                            "^^ 

576 

C 
C 

ELEVATION  SUBROUTINE  ^ 

DIMENSION  STAIlboTrELEVI 100) .XLVCIIOJ.PCSTAIIOI f PTSTA(iO) 

c 

READ  FORMATS 

577 

900 

FORMAT(2F10.2)                         _.    -.  .-  .  -      --   . 

578 

c 
c 

901 

F0RMAT(F10.2) 

WRITE  FORMATS                        ...     ^ ...^  -   - 

579 

953 

F0RMAT(33X,15HBEGINNING  POINT/) 

580 

954 

F0RMAT(31X,9HSTATI0N  =, Fll . 2/31X, IIHELEVAT ION  =,F9.2//) 

581 

955 

F0RMAT(36X,9HEr4D  POINT/) 

582 

960 

1 

F0RMAT136X,5HCU^VE,I3//30X,I4HCURVE  LENGTH  =, F9. 2/31X, 12HPC  STATIO 
IN  =,F10.2/31X,12HPT  STATION  =,F10.2//) 

583 

JN=NPVI*2                             .  . 

.  58^ 

NP=NPVI*1 

585 

KK=NNN*l 

586 

KKK=NN\»2»NPVI 

c 

END  POINTS  MUST  BE  BEYOND  STATIONED  POINTS 

c 

♦ 

READ  STATIONS  AND  ELEVATIONS  OF  PVI'S  INCLUDING  END  POINTS 

587 

00  10  K=KK,KKK                                    ... 

588 

10 

READ  (1,900)  STA(K),ELEV(K) 

c 

THE  END  PVI'S  ARE  ASSIGNED  A  CUf^VE  LENGTH  OF  ZERO  SO  THAT  THE  FIRS 

c 

AND  THE  LAST  PVC  WILL  NOT  BE  ON  THE  HORIZONTAL  ALIGNMENT  

589 

XLVCI l)=0. 

590 

XLVC(JN)=0. 

591 

!F(NPVI)  30,30,15                     .              _.   ^       

c 

♦ 

ASSIGN  LENGTHS  OF  VERTICAL  CURVES  AT  EACH  PVI 

592 

15 

DO  20  J=2,NP 

593 

20 

READ  (1,901)  XLVC(J) 

c 

COMPUTE  STATIONS  OF  PC'S  AND  PT'S 

594 

30 

00  40  K=KK,KKK 

595 

C2=XLVC(K-NNN)/2.              

596 

PCSTA(K-NNN)=STA(K)-C2 

597 

PTSTA{K-NNN)=STA(K)+C2 

598 

40 

CONTINUE                                '                       

c 

WRITE  STATION  AND  ELEVATION  OF  BEGINNING  POINT 

599 

WRITE  (3,953) 

600 

WRITE  (3,954)  S  TA  (  KK)  ,  ELEV  (  KK)          .^  

601 

IF(NPVI)  46,46,42 

c 

WRITE  VERTICAL  CURVE  GEOMETRY 

602 

42 

DO  45  1  =  2, NP                        .    _               ._^ 

603 

K=I-1 

604 

45 

WRITE  (3,960)  K , XLVC( I ) , PCSTA ( I ) , PTSTA ( I ) 

c 

WRITE  STATION  AND  ELEVATION  OF  END  POINT       _  

605 

46 

WRITE  (3,955) 

606 

WRITE  (3,954)  S T A ( KKK ) , ELEV (KKK ) 

\ 

c 

COMPUTE  ELEVATIONS  AT  EACH  POINT                .     

607 

DO  100  1=1, NNN 

6C8 

DO  60  K=KK,KKK 

609 

IF(PCSTA(K-NNN)-STA(I ))50,70,70 

610 

50 

IF(PTSrA(K-NNN)-STA(I )) 60, 60, 80 

611 

60 

CONTINUE 

c 

POINT  ON  TANGENT 

612 

70 

ELEV(I)=(STA(I)-STA(K-1))»(ELEV(K)-ELEV(K-1))/(STA(K)-STA(K-1)) 

613 

ELEV( I)=ELEV(I) *ELEV{K-1) 

614 

GO  TO  100 

c 

POINT  CN  VERTICAL  CURVE 

615 

80 

ELEV( I)  =  (STA( I)-STA(K-l) )*( ELEV ( K )-ELEV ( K-1 ))/( STA( K ) -STACK- 1 J ) 

616 

ELEV(I)=ELEV(I)*ELEV(K-1) 

i 

617 

618 
619 
620 
621 
622 
623 
624 


100 


G2=  {  ELE  V  (  KU  )  -EL  EV  (  K  >  >  /  {  STA  ;  K  M  )  -ST  A  {  K  )  ) 
G1=(ELEV(K)-ELEV(K-I) ) / ( STA ( K) -STA ( K-l ) ) 
Y=(G2-Gl)*(SrA(  I  »-PCSTA(K-NNN)  )*(SrA(I  j-PCSTACK-NNiNU 
Y=Y/(XLVC(K-NNiN)*2.  ) 
ELEV(  I)  =  ELEV(n+Y 

C0^1TI^UE  ._ 

RETURN 
END 


5  7 


630  SUBROUTINE    IrNTERP     (  Ai  ,  AEDGEi  A2,  BltB2,C2) 

631  C2=B1+(B2-BI)*(AEDGE-Al)/(A2-Al) 

632 RETURN  _  ^      

633  END 


58 


634  SUB'-ICUTINE  Df'xAW  {  L,  N»}!.  V,  NSK  IP  } 

635  DIMENSION  NSKIP ( IS) ,11  ( 15 , 300) , V ( 1 5, 300 ) , TNPIC ( 15 , 300) 
C      ISKIP  IS  TWO  TIMES  THE  INTERVAL  NUMBER  FOR  CROSS  LINES. 

636  LSAVE=L 

637  ISKIP=NSKIP(17)*2 
638 WRITE  (3,505)  _- 

F0RMAT(7(/) ,16H  DRAW  SUBROUTINE) 
HMAX=10,0 

HMIN=-10.0 

VMAX=5.0 
VMIN=-5.0 

PROGRAM  DRAWS  LINES  BACK  AND  FORTH  

DO  100  1=1, L, 2 
IH1  =  0 

IV1  =  0      _  _..   . 

NFRST=0 
IHV3=0 

NN=NSKIP{ I )  

IF(NN.EQ.O)  GO  TO  46 
DO  7  K=1,NN 

I NP I C  (  I  ,  K  )  =1 : 

NN=NN+1 

DRAWS  LINES  1,3,5,7  FROM  1  TO  N 

DO  98  J=NN,N      

NEXT  POINT  H2,  V2 
H2=H( I, J) 

V2=V { I , J )   

!NPIC{I,J)=0 

IH2  =  0 

HMI N  CHECK         

659  IF(H2-HMIN)1,2,2 

660  1  IH2=1 

661 INPIC{I,J)  =  1    

*  662  IF(NFRST)2,18,2 
,  663        2  IF( IH1-1)3,3,4 

664 3  IF(IH2-IH1)5,4,6       ..       

665  5  CALL  INTERP ( HI, HMIN,H2, VI , V2, VI ) 

666  H1=HMIN 

^  667  IF(V1-VMIN)31,9,8 

*  668        8  IF(VMAX-V1)47,9,9 
,   669       31  IVl=l 
__.670 GO  TO  10     

671  47  IV1=2 

672  GO  TO  10 
.„673 „_   9  IVI  =  0 

674       10  IP1=1 
»  675          GO  TO  18 
676    „    6  H3=H2         _  _   

677  V3=V2 

678  IHV3=1 

679   CALL  INTERP{H1,HMIN,H2,VI,V2,V2) 

680  H2=HMIN 

*  681          IH3=IH2 
_-682 IH2  =  0      .^  

683  GO  TO  18 

C      HMAX  CHECK 

^—684 4  IF(HMAX-H2)11,12,12 

685       11  IH2=2 

*  686         INPICCI,J)=1 
--687 ...  IFINFRST)12,18,  12     


639 

505 

640 

'  641 

642 

643 

C 

644 

.  645 

646 

647 

648 

649 

650 

,  651 

__652^ 

7 

653 

46 

C 

_654  _ 

C 

V  655 

__656  _ 

657 

658 

C 

59 


688 

12 

-  689 

13 

690 

691 

692 

15 

693  

J»8 

694 

>  695 

49 

696 

697 

17 

698 

16 

699  

700 

T4 

*  701 

702 

703 

704 

705 

706 

►  707 

18 

C 

708 

709 

19 

710 

711 

»  712 

20 

713  „_ 

21 

714 

41 

715 

23 

716 

717 

718 

.719 

42 

720 

24 

'  721 

722 

723 

724 

»  725 

726 

C 

^727   _ 

22 

728 

25 

729 

i_730 

731 

26 

>  732 

43 

.733 

27 

734 

735 

i     736 

737 

44 

^  738 

28 

739 

740 

741 

»  742 

743 

'   744 

„J45 

29 

IF(IH2-IH1>13,18,14 

CALL  INTeRP(Hl,HMA>C,H2,Vl,V2,Vn 

H1  =  HMAX . 

IF(V1-VMIN)48, 17,15 
IF{VMAX-\/l)49,17,l7 

IV1  =  1 

GO  TO  16 
IVl  =  2 

GO  TO  16 : 

IV1  =  0 
IP1=1 

GO  TO  IB 

H3  =  H2 

V3  =  V2 

IHV3=1 

CALL  INTERP(H1,HMAX,H2,V1,V2,V2) 

H2=HMAX 

IH3=IH2     

IH2  =  0 
IV2  =  0 

VMIN  CHECK  '. 

IF(V2-VMIN)19,20,20 
IV2=1 

INPIC(I,J)=1         

IF(NFRST)20,37,20 

IF(IVl-l)2l, 21,22 

IF(IV2-IV1)41,22,42 

IF(  IH2)37,23,37 

CALL  INTERP(V1,VMIN,V2,H1,H2,H1) 

V1=VHIN 
IPl^l 
GO  TO  32 

IF(IH2)37,24,37      .  .  _. 

H3  =  H2 
V3--V2 

IHV3=l  .  - 

IH3=IH2 

CALL  INTERP(V1,VMIN,V2,HI,H2,H2) 

V2=VMIN  

GO  TO  32 
VMAX  CHECK 

IF(VMAX-V2)25,26,26  _ „ 

IV2  =  2 

INPIC( I,J)=1 

IF(NFRST)26»37,26      .  _ 

IF(IV2-IV1)43,29,44 

IF(IH2)37,27,37 

CALL  INTERP(V1,VMAX,V2,H1,H2,H1) 

V1=VMAX 

IP1  =  1 

GO  TO  32  _._   ... 

IFl  IH2)37,28,37 
H3=H2 

V3=V2  .  . 

IHV3=1 
IH3=IH2 

CALL  INTERPCVl, VMAX,V2,H1,H2,H2) 
V2=VMAX 
GO  TO  32 
I  IF( IH2)37,30,37     


60 

746  30  ir-nv2)37,32,37 

C      ORIGINAL  POINT  OUT  OF  PICTURE 

747  32  IF{NFRST)34,33,34                               . 

748  33  CALL  PLOT (-V2,+H2, 3 ) 

749  GO  TO  37 

750  34  IF{  I  PI)  36, 36, 35 

751  35  CALL  PLOT (-V1 ,+Hl,3 ) 
»  752       36  CALL  PLOT ( -V2,  +  H2 ,2  ) 

753       37  IF( IHV3)38,38,39 

._ 

~754'"     38  Hl  =  H2 

755  V1=V2 

756  GO  TO  40                                  _  .   

757       39  Hl=H3 
*  758          V1=V3 
159          IH2=IH3 

760  40  NFRST^l 

761  IH1=IH2 

762  .   ___    IV1=IV2     .           .   „     -   ..   

763  IHV3=0 
».  764          IP1  =  0 

765       98  CONTINUE 

: 

C      DRAWS  LINES  2,4,6,8  FROM  N  TO  1 

766  NFRST=0 

767  11=1+1 

768          IFdl.GT.L)  GO  TO  100 
r  769          IH1=0 
770          IV1=0 

771  NN=NSKIP(I+1) 

772  J  =  Nn 

773  IF(NN.EO.O)  GO  TO  97 

774  DO  197  K=1,NN 

775  197  INPIC( I+1,K)=1 

776  97  NN=NN+1 



777          DO  99  K=NN,N 
r     778          J=J-1 
.  .  C      NEXT  POINT  H2,  V2 

779  H2=H(I+1,J) 

780  V2=V(I+1,J) 
^.781          IiNPICII  +  l,J)=0 

782  IH2=0 

C      HMIN  CHECK 

783  IF(H2-HMIN)51,52,52 

784  51  IH2=1 

785  INPIC{I+1,J)=1 
^..786  _ ^     IF(NFRST)52,68,52 



787       52  IF(IH1-1)53,53,54 
»   788       53  IF(IH2-IH1)55,54,56 
._  789 55  CALL  INTERP  (  HI,  HMI  N,H2,  VI ,  V2,  VI ) 

790  Hl=HMIN 

791  IF(V1-VMIN)50,59,58 
»..792.  58  IF(VMAX-V1)81,59,59 

793       50  IV1=1 
^   794          GO  TO  60 
_  795 81  IV1  =  2 



796  GO  TO  60 

797  59  IV1=0 
k___798 60  IP  1=1 

799          GO  TO  68 
'  800       56  H3=H2 
801  V3=V2 

802 

IHV3=1 

803 

CALL  INTERP(Hl,HMIN,li2, 

iVl, 

»V2, 

,V2) 

9 
1 

804 
805 
806 



H2  =  HiMIN 
IH3=IH2 
IH2  =  0 

, ,-.  ... 

_807^_ 

C 

GO  TO  68 
UMAX  CHECK 





808 

54 

IF(HMAX-H2)61,62,62 

» 

809 
810 
811 

61 

IH2=2 

INPIC(I+l,J)=l 

IF(MFRST)62,68,62 





812 

62 

IF(IH2-IH1)63,68,64 
CALL  INT£RP{H1,HMAX,H2, 

rVL 

tV2, 

813 

63 

,V1) 

814 

H1=HMAX 

^ 

815 
816 
817 

65 
96 

IF(V1-VMIN)96,67,65 

IFtVMAX-Vl)118,67,67 

IV1=1 

-  ■  ■ 

818_ 
819 

820 

118 

GO  TO  66 
IV1  =  2 
GO  TO  66 

— 



821 
822 
823 

67 
66 

IV1  =  0 
IP1  =  1 
GO  TO  68 







_824_ 
825 
826 

64 

H3=H2 
V3=V2 
IHV3=1 

> 

827 

828 
829 



CALL  INTERP(H1,HHAX,H2, 

H2=HMAX 

IH3=IH2 

»V1 

,V2, 

rV2) 

830 

IH2  =  0 

831 

68 
C 

IV2=0 
VMIN  CHECK 

832 

833 

69 

IF(V2-VMIN)69,70,70 
IV2=1 





4 

834 

INPICII+1,J)=1 

835 

836 
837 
838 

'  70 

71 
91 

IF(NFRST)70,87,70 
IF( iVl-l)71,71,72 
IF( IV2-IV1)91,72,92 
1F(IH2)87,73,87 

._ 

V 

839 
840 

73 

CALL  INTeRP(VUVMIN,V2, 
V1=VMIN 

rHl 

iH2, 

,H1) 

-— 

841 
842 
843 

92 

IPI=l 

GO  TO  82 
IF(IH2)87,74,87 

-- 

... 

844 

74 

H3=H2 

t 

845 

V3=V2 

1 

846 

IHV3=1 

• 

9 

847 

IH3=IH2 

848 

CALL  INTERP(V1,7MIN,V2, 

rHl. 

»H2, 

rH2) 

849 

V2=VMIN 

850 

GO  TO  82 

^ 

C 

VMAX  CHECK 

\ 

851 

72 

IF(VHAX-V2)75,76,76 

852 

75 

IV2  =  2 

853 

INPIC(I+1,J)=1 

854 

IF(NFRST)76,87,76 

855  _ 

76 

IF[IV2-IVI)93,79,94 

> 

856 

93 

IF(IH2)87,77,87 

V 

857 

77 

CALL  INTERP(VI,VMAX,V2, 

rHl, 

»H2, 

,H1) 

858. 

Vl=VMAX 

61 


62 


859 
860 
861 
862 
863 
864 
863 
866 
867 
868 
869 
,870 

871 
872 
873 
874 
875 
876 
877 
878 
879 
880 
881 
882 
883 
884 
885 
886 
887 
688 
889 
890 
891 
892 
893 
894 
895 

896 
897 
898 
899 
900 
901 
902 
903 
904 
905 
906 
907 
9C8 
909 
910 
911 
912 
913 
914 
915 


94 
78 


79 

80 

82 
83 

84 

85 
86 
87 
88 


89 


90 


99 

100 

148 

501 

502 

147 
140 


143 

144 

145 
146 

141 


IP1  =  1 

GO  TO  82 

1F(IH2)87,78,87._  .__    _  __  .     ^   _  ^ 

H3=H2 

V3:=V2 

IHV3=l  .  .^ _      

IH3=IH2 

CALL  INTERP(V1,VMAX,V2,H1,H2,H2) 

V2=VMAX     ^_  ^         __ 

GO  TO  82 

IF(IH2)87,80,87 

IF( IV2)37,8?,87  _ 

ORIGINAL  POINT  OUT  OF  PICTURE 

IF(NFRST)84,83, 84 

CALL  PLOT(-V2,+H2,3)    _  

GO  TO  87 
IFnPl)86,86,85 

CALL  PL0T(-V1,4-H1,3)      , 

CALL  PL0T(-V2,+H2,2) 
IF( IHV3)88,88,89 

H1  =  H2 __•    _   

V1  =  V2 

GO  TO  90 

Hl  =  H3    . 

vi=y3 

IH2=IH3 

NFRST=l  .  _  _         

IH1=IH2 
IV1=IV2 

IHV3  =  0     _        _  „   _    

IP1  =  0 
CONTINUE 

CONTINUE  _      

IF(NSKIP(18) )  147,147,148 

WRITE(3t50l) 

FQRMAr(lHI,6(I3H   I   J   INPIC,8X)//) 

WRITE  (3,502)  ( ( I , J, I NPIC ( I , J ) , J= 1,N) , 1=1 ,L) 

F0RMAT(6(lX,2I3,r5,9X)  ) 


IFtNS 

L  =  4 

00  14 

GO  TO 

JJ=1 

KK=2 

GO  TO 

JJ=3 

KK^-1 

GO  TO 

JJ  =  2 

KK=3 

NN  =  NS 

DO  14 

INPIC 

HIKK, 

V(KK, 

NSKIP 

NSKIP 

NSKIP 

DRAWS 


KIP(16))  140,140,142 

1  J=l,3 
(143,144, 145), J 


146 


146 


KIP{JJ)+1 

1  I=NN,N 

(KK,I)  =  INPIC{JJ,  I) 

I)=H(  jj,n 

I)=V( JJ,I) 
(2)=NSKIP{ I) 
(l)=NSKIP{3) 
(3)=NSKIP(2) 
CROSS  LINES  AT  EVERY 


STATION  BACK  AND  FORTH 


63 


916 

142 

917 

_  _504 

918 

919 

920 

921 

102 

922 

101 

923 

924 

925 

300 

926 

927 

928 

929 

930 

931 

932 

933 

934 

935 

936 

937 

938 

939 

9^0 

941 

942 

943 

944 

945 

946 

947 

948 

949 

950 

951 

952 

953 

954 

955 

956 

301 

C 

957 

302 

958 

959 

960 

961 

c 

962 

963 

964 

965 

103 

966 

105 

967 

968 

969 

970 

971 

106 

972 

WRITE  (3,504) 

FORMATCM/)  ,12H  CROSS  LINES)      _____   _ 

HSKIP=NSKIP(1) 
DO  101  K=2,L 

IF{NSKIP(K)-MSKIP)l02,101f  101  __; 

MSKIP=NSKIP(K) 
CO^NTINUE 

MSKIP  =  MSKIP4-1  _     

IF(NSKIP(16)-4)  302,300,302 

L^14 

DO  301  1=1, N  _   _ 

H(14,I)=H(12,n 

V(14,I)=V(12,I) 

H(13,I)=H(11,I)  __^    „ 

V(13, I)=V(11,I) 

H(12, I)=H(10,I) 

V{12,I)=V(10,I)  „   _     _  ___ 

H{ll,I)=Hl9,n 

V(11,I)=V(9,I) 

H{10,I)=H(9,I)     • 

V(10,I)=V{9,I) 
H(9,I)=H(8,n 

v(9,r)=v(8,n    _ 

H(8,I  )  =  H(7,I  ) 
V{8,I)=V(7,I) 
H(7, I)=H(6,I) 
V(7,I  )  =  V(6,I) 
H(6,I  )  =  H(5,I) 

V(6,I  )  =  V(5,  I  )         _  

H(5,I)=H(4,I) 
V(5,I  )  =  V{4,  I  ) 
INPIC(14,I  )  =  INPIC(12,I)      _         _____      ___   __ 

INPIC{13,I )=INP[C(11,I) 
INPIC(12, I )=INPIC(10,I) 

INPIC(I1,I)  =  INPIC(9,I  )       __   _    _ 

IMPICdCI  )  =  INPIC(9,I)  " 

INPIC{9,I)=INPIC{8,I) 

INPIC(8,I)=INPIC(7,I)        _ 

INPIC(7,I)=INPIC(6,I) 

INPIC(6,I)=INPIC(5,I) 

INPIC15, I)=INPIC(4,I) 

STARTS  CROSS  LINES  AT  BOTTOM  OF  PICTURE,  PROGRESSES  TO  TOP 

DO  200  JJ=MSKIP,N, ISKIP 

J=JJ  _     _  _ 

NFRST=0 

ID=1  .     . 

I=L+1  

DRAWS  MSKIP,  MSKIP+2,ETC,  LINES  RIGHT  TO  LEFT 

00  198  K=l,L,2 

1  =  1-1  _  _    ^        _   _ 

IF(NFRST) 104,103,104        '  " 

IF(INPIC( I, J) )106,105,106 

SAVEV=-V( I, J)  _      

SAVEH=H(I,J)  '    " '    '  "        ' 

KH=I 

KV=J  

10=0  '     "^       "" ■  ■  " 

NFRST=1 
GO  TO  112 


973 

lO^i 

974 

107 

>■    975 

^109 

976 

977 

978 

979 

980 

•  981 

108 

,  982 

110 

983 

984 

.  Ill 

985 

•  986 

112 

'  987 

»  988 

114 

989 

116 

.  990  . 

991 

117 

"  992 

198 

L_993 

♦  994 

995 

_._996_ . 

997 

* 

C 

'  998^ 

.  999 

1000 

1001 

153 

1002 

155 

1003 

.  1004 

.  1005 

1006 

'1007 

156 

10C8 

1009 

154 

.  1010 

157 

■  1011 

159 

1012 

1013 

1014 

1015 

.  1016 

'  1017 

158 

1018 

160 

1019 

1020 

161 

1021 

•  1022 

162 

1023 

.  1024 

164 

1025 

166 

1026 

1027 

167 

'  1028 

199 

1029 

200 

,  1030 

1031 

IF( 10)107,103,10? 

IF(  INPICd,  J)  )112,109,112 

SAVEv=-v(  [,j)  .. ; 

SAVcH=H{I,JJ 
KH=I 

Kv=j  _  _     : 

ID=0 

GO  TO  112 

IF(INPIC(I,J) )110,111,110         __   _ 

I  0=1 

GO  TO  112 

CALL  PLOT  (SAVEV,SAVEH,3)         _._    ; 

CALL  PLOT(-V{I, J),+H{I,J),2) 
1  =  1-1 

IF{  10)198,114,198  __ _ 

IFdNPICt  I,  J)  )117,116,I17 
CALL  PLOT  (SAVEV,SAVEH,3) 

CALL  PLOT(-V( I, J),+H(I,J),2)     

10=1 
CONTINUE 

J  =  JJ  +  NSKIP(17)        J 

IF(J.GT,N)  GO  TO  200 
NFRST=0 

10=1   .„ _„ 

1=0 

DRAWS  MSKIP+1,  MSKIP+3,  ETC,  LINES  LEFT  TO  RIGHT 

DO  199  K=1,L,2 _  ^ 

1  =  1  +  1 
IF(NFRST)154,153,154 

IF( INPIC( I,J) )156,155,156  _    _  

SAVEV=-V{ I, J) 
SAVEH=H(I, J) 

KH=I     ._ 

KV=J 
ID=0 

NFRST  =  1  _ 

GO    TO    162 
IF(ID)157,158,157 

IF(INPIC{I,J) )162,159,162  

SAVEV=-V( I, J) 
SAVEH=H(I, J) 

KH=I 

KV=J 
10=0 

GO   TO    162  _  _      _ 

IFdNPICd.J)  )160,161,160 
ID=l 

GO    TO    162 

CALL  PLOT  (SAVEV,SAVEH,3) 
CALL  PLOT(-V(I,  J),+Hd,J),2) 

1  =  1  +  1  ^ 

IFdOigg,  164,199 
IFdNPICd.J)  )167,166,167 

CALL    PLOT    (SAVEV,SAVEH,3)  _      ^       _      _ 

CALL    PLOT(-Vd,  J),+Hd,J),2) 
ID=1 

CONTINUE  .       ^  _ 

CONTINUE  '"    "  '"    • "  ■ 

L=LSAVE 

WRITE    (3,506)    HSKIP 


b'j 


1032  506  FORMAT  (//6X, 'FIRST  POINT  ***',I';} 

1033  RETURN 

1034  END 
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ABSTRACT 


The  purpose  of  this  study  was  to  investigate  the  impor- 
tance of  the  spiral  curve  in  the  visual  appearance  of  a 
horizontal  curve.   This  was  to  be  accomplished  by  simulating 
actual  and  theoretical  conditions  of  roadway  geometry  and 
rating  the  visual  appearance  of  each  location.   The  simulation 
involved  converting  three-dimensional  coordinates  into  two 
dimensional  perspective  coordinates  and  plotting  these 
coordinates  thus  giving  a  perspective  drawing  of  the  roadway. 

Many  different  coml^inations  of  sight  distance,  display 
angle  and  roadv/ay  geometry  were  simulated  and  rated  in  an 
attempt  to  determine  the  factors  which  affect  the  visual 
appearance  of  a  roadway.   It  was  found  that  increasing  sight 
distance  caused  the  appearance  of  a  curve  to  become  less 
acceptable.   The  display  angle  was  found  to  be  proportional  to 
the  appearance  of  the  curve,  i.e.  an  increase  in  the  display 
angle  resulted  in  an  improved  visual  appearance  of  the  curve. 
The  geometry  of  the  curve,  spiral  length,  likewise  affected 
the  visual  acceptability  of  the  curve.   The  longer  the  length 
of  spiral  used,  the  more  visually  acceptable  it  became. 

A  preliminary^  investigation  was  undertaken  to  determine 
the  feasibility  of  using  the  rate  of  change  of  slope  of  a 
curve  in  the  perspective  picture  plane  as  a  means  of  determin- 
ing the  visual  appearance  of  a  curve  v/ithout  drawing  a 
perspective  view.   This  investigation  indicated  that  this  approach 
does  warrant  further  study. 


